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A recurring motif in the organization of biological tissues are networks of long, fibrillar protein 
strands effectively confined to cylindrical surfaces. Often, the fibers in such curved, quasi-2D 
geometries adopt a characteristic order: the fibers wrap around the central axis at an angle which 
varies with radius and, in several cases, is strongly bimodally distributed. In this Letter, we 
investigate the general problem of a 2D crosslinked network of semiflexible fibers confined to a 
cylindrical substrate, and demonstrate that in such systems the trade-off between bending and 
stretching energies, very generically, gives rise to cross-hatched order. We discuss its general 
dependency on the radius of the confining cylinder, and present an intuitive model that illustrates 
the basic physical principle of curvature-induced order. Our findings shed new light on the potential 
origin of some curiously universal fiber orientational distributions in tissue biology, and suggests 
novel ways in which synthetic polymeric soft materials may be instructed or programmed to exhibit 
preselected macromolecular ordering. 

PACS numbers: 87.15.A-, 82.35.Pq, 82.35.Gh 


In the arterial wall, quasi-2D layers of meshed colla¬ 
gen fibers alternate with layers of smooth muscle cells 
to envelop and strengthen the cylindrical lumen of the 
blood vessel pQ. Likewise, in the annulus hbrosus of the 
intervertebral disc concentric, thin lamellae of collagen 
cylindrically surround the soft nucleus pulposus [2H3] , an 
ordering that is proposed to find its origin in the collagen 
sheet that wraps around the notochord, the embryonal 
precursor of the spine [5]. Elsewhere in the body, in corti¬ 
cal bone thin lamellae of mineralized collagen are wrapped 
cylindrically around each Haversian canal in a structure 
called the osteon iMj. In each of these, physiologically 
unrelated, settings a universal organization is employed: 
thin networks of long, fibrillar protein strands tightly 
enveloping cylindrical domains. Remarkaby, in each of 
these cases, fibers in such curved, quasi-2D geometries 
tend to become ordered: the biopolymers wrap around 
the central axis at an angle which varies with radius and, 
in several cases, is strongly bimodally distributed, leading 
to a cross-hatched appearance. Several examples of such 
order are depicted in Fig. |Tj 

Order in filamentous biomaterials may be provoked in 
several manners. In |9j, it is reported strain may induce 
alignment. Part of this alignment is imprinted irreversibly. 
In [TUIITT] , it is demonstrated that strain, likewise, may 
affect the enzymatic degredation of collagen in living tis¬ 
sues. Strain-dependent fiber deposition and degradation 
- mimicking cellular or enzymatic activity - is included 
in several constitutive models HM2, and the interac¬ 
tions between cylinders (carbon nanotubes and histones, 
in particular) and polymers (either single polymers, or 
dilute solutions without linking) have been adressed in 
several previous works, reporting helical wrapping con¬ 
figurations - but in these cases the polymers were not 
fully confined to the surface [MU! but rather interacted 



FIG. 1: (a) Cross-hatched order of collagen in the intima of the 
human thoracic aorta JT[. (b) The orientational distribution 
function for the angle <j> with the cylindrical axis of the fibers 
in image (a), (c) Cross-hatched ordering of collagen in the 
midsagittal section of the annulus fibrosus [I]. 


with it through a weak attractive potential. Though these 
models all display some degree of ordering, they fail to 
provide a mechanism for this ordering in the absence 
of active orienting entities such as cells - the manner in 
which it appears to occur in many biological systems. We 
do not claim that in all of the biological instances we 
have mentioned a single, physical mechanism explains all 
of the observed structure (indeed, in several cases the 
cross-hatched order is absent, or replaced by more helical 
patterns). What is, however, clear is that the interaction 
between curved surfaces and semiflexible polymers, in 
general, gives rise to nontrivial ordering effects that are 
manifested in several, biologically relevant, settings. 

In order to understand and, perhaps, learn how to 
control this ordering we investigate a simple model of 
a polymer meshwork fully confined to the surface of a 
cylinder. 

Model. Our polymer chains are confined to the surface 
of a cylinder. We shall denote by (u, v) the two coor¬ 
dinates on the rolled-out cylinder with 0 < u < L the 
axial coordinate, and 0 < v < 2 tt the circumferential 
coordinate; these coordinates are mapped in 3D onto 
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the surface of a cylinder of radius R and length L lying 
along the y-axis: r(u,v) = {x(u, v), y(u, v), z(u, i>)} = 
{Rcos(v/R),u,Rsin(v/R)}. The semiflexible polymers 
themselves are modeled as discrete, extensible persistent 
chains. That is, chains of total, unstretched contour 
length ft are represented by discrete segments labeled 
by * = !,■■■ ,N, each such segment having a differ¬ 
ent rest length (ft. Segment i follows bond 2-vector 
ft = {uj+i — tq, |_i — Vi} on the rolled-out cylinder, and 
bond 3-vector /3j = r(tt,; + i, i>j + i) — r(rq, tft in 3D. The 2- 
and 3-lengths of segments are denoted by ii = (bi ■ bi) 1 / 2 
and Xi = {fa ■ Pi) 1 / 2 , respectively. All segments possess 
the same extensional modulus Y, and correspondingly 
have different spring constants ft = Y/ft^. Stretching 
of each segment is modeled as a simple Hookean spring, 
with a stretching energy 'H stle tch,i = 5 ^ 7 (^i - ft,*) 2 - For 
the bending energy of the chain, there are two contribu¬ 
tions that one must consider. The first is the bending 
between segments, proportional to the square of the angle 
between two adjacent segments. The second, however, 
arises uniquely on curved substrates and measures the 
bending of each segment individually. First, a note on the 
intersegment bending term. This term, proportional to 
the cosine of the angle between subsequent segments, may 
be written in terms of the angle ft between adjacent bond 
vectors /ft and ft +i as 7ft nte rseg,j = a,+a 1+ i C 1 “ cos ft). 
However, we will consider here networks fully confined to 
the cylindrical surface, and will not allow for excursions 
or shortcuts through the interior. That is, we will assume 
that li = A,;, as the mapping onto the cylinder surface is 
isometric. A consequence of this is that each individual 
segment is bent to conform to the surface - the additional 
energy this imparts is treated later on. We may now 
write the angle bending term, instead, as a function of 
the angle ft between segments in the (it, v) plane. Due to 
the isometric mapping the in-plane angle 9 at which two 
bent segments join remains the same as the 3D angle H 
(this is no longer true for general curved surfaces). Thus, 
our angle bending potential is chosen as 

Wi n terseg,{i,i+1} = ~n j ~ 0 (1 — COS ft) ~ — - - ft , 

+ ft- 1-1 ft + 1 

(1) 

where in the final equation we assume small in-plane 
angles between the fibers. As stated, the energy of the 
chain picks up an additional segment contribution, which 
measures the bending of segments to conform to the sur¬ 
face. This we compute by first constructing the geodesic 
7 j(s), connecting the points (iq,i>j) and (itj + i,Uj+i), with 
0 < s < ii its natural (arclength) parameter. The 
geodesic is a simple linear form in s in this case, and 
the curvature contribution is then straighforwardly (and 
analytically) evaluated to be 


Wsegbend,?' — ^ 


r 

d 2 ft 

lo 

ds 2 


dS= 2 W^i), ( 2 ) 


with pi the angle between the segment and the y-axis (i.e., 
the axial direction of the cylindrical substrate): cos (ft) = 

ft-ft 

The appearance of this angle should come as no surprise 
- a segment aligned to the long axis of the cylinder (</> = 0) 
remains unbent, whereas one that is circumferentially 
oriented (</> = ±7r/2) is maximally bent, when wrapped 
around the cylinder. This is the origin of the effects 
we discuss here: the curved substrate acts, locally, as 
an external field in the in-directions. Obviously, the 
term vansishes as R~ 2 when the radius of the cylinder is 
increased. While one might be tempted to expand the 
expression for "H seg bend around <f> = 0, this is not justified 
as the random arrangement of chains in a crosslinked 
network consists of segments oriented at all angles with 
respect to the cylinder’s long axis. Thus, our complete 
energy functional is 



In what follows, we attempt to characterize the typi¬ 
cal spatial arrangement of the fibers in networks that 
minimize this energy functional. Before we move on to 
crosslinked networks, it is instructive to first consider the 
highly simplified model of a single ring polymer wrapping 
around a cylinder. 

Ring Polymer. For a simple, continuous ring polymer 
wrapped around a cylinder of radius R, the problem and 
the typical nature of solutions becomes clear. We will 
assume, again, that the polymer has a rest length ft, a 
spring constant Y/ft and a bending modulus n. The trade¬ 
off for this ring polymer, in terms of energy contributions, 
is obvious: circumferential positioning allows the polymer 
minimal length, but at the expense of maximal curvature 
(since the radius R is the only curvature length scale in 
the problem). Slanted arrangements, where the polymer 
traces out an ellipse that forms the intersection between 
the cylinder and a plane at some angle, require less cur¬ 
vature energy but more stretching. We call <fr the angle of 
the ring polymer with respect to the cylinder’s long axis 
(i.e., <j> = ±7t/ 2 corresponds to circumferential placement, 
and (j> = 0 to infinitely stretched out axial placement), 
and call (ft the corresponding elliptic conformation. Since 
the polymer remains planar, (ft may be parametrically 
represented in 2D as <ft(f) = i?{sin _ (</!>) cos(f), sin(f)}; 
0 < t < 2n. We will denote by /ft the arc length corre¬ 
sponding to one trip around the cylinder at angle 4> (a 
closed form expression in terms of elliptic integrals exists 
for /ft, but this is not particularly instructive). The opti¬ 
mal angle </>, given Y and k, may be determined now by 
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FIG. 2: (a) The ring model: a semiflexible ring polymer (blue) 
is wrapped around a cylinder of radius R, at an angle <j> with 
respect to the cylinder’s long axis. A competition between 
bending (favoring axial alignment) and stretching (favoring 
circumferential alignment) results in a bifurcation, (b) The 
optimal angle </> m i n as a function of the radius for a polymer 
of length Lo = 2. 


minimizing over all <f> the energy 


^ring(</>) 


2^ _£o)2 + 



d s. 


(4) 


Evaluating this expression, and minimizing it as a func¬ 
tion of <j> yields an interesting result, depicted in Fig. [2] 
For large radii of the cylinder, the bending terms become 
unimportant and only the stretching remains. The stretch¬ 
ing terms simply seeks to minimize the contour length, 
and indeed the energy is minimized at <f> = 7t/ 2; the ring 
encircles the cylinder orthogonal to its long axis. If the 
radius is lowered, however, beyond a critical radius R*, 
the bending becomes increasingly important and the cir¬ 
cumferential ordering gives way to two related slanted ori¬ 
entations as the new minima. The critical radius R* may 
be expressed in dimensionless form: R */£ 0 = ^(Y^q/k) 
with ip(x) a universal, though rather unwieldy, function of 
x. The upshot of this simple ring model is, that for closed 
filaments confined to the surface of a cylinder, below a crit¬ 
ical radius the circumferential ordering is superseded as 
the minimal energy configuration by two slanted orienta¬ 
tions, whose angles <fi are symmetric around tt/2. Clearly, 
this reminds us of the crossed order that is seen in the bi¬ 
ological systems we recalled in the introduction. We now 
turn to simulations of fully crosslinked networks of poly¬ 
mers to verify to what extent this behavior is manifested 
in more complex settings. 

Network simulations. To simulate the behavior of many 
polymers, crosslinked to each other and confined to the 
surface of a cylinder, we adopt the mikado procedure: 
an initial network of straight fibers is prepared by ran¬ 
domly distributing them on the rolled-out cylinder (i.e., 
in the (it, u)-plane). The chains interact with each other 
only at crosslinking points, there is no excluded volume 
interaction. Dangling ends are removed. In this initial 
network, the energy contributions from stretching "Hstretch 
and bending "Hmterseg are zero by construction. The seg¬ 
ment bending term ’H seg bend, however, is not. While any 
such mikado arrangement is a minimizer of the energy 
functional 'H for R = oo, this is not true for finite radii. 


Therefore, we allow the initial network to relax to a mini¬ 
mal energy state, in which some of the segment bending 
energy is alleviated by reorienting the segments. The 
elastic energy lost in doing so is, of course, partly traded 
for stretching and intersegment bending energies. To find 
the energy minimum, we employ a multistep non-linear 
conjugate gradient scheme m- 

The result of such a minimization is shown in Fig. [3] 
Figs 3(al), 3(a2) and 3(a3) display data from the initial, 
disordered configuration, where Hsegbend is the only non¬ 
zero energy component. Figs. 3(bl), 3(b2) and 3(b3) 
present the network as it appears after the minimization 
procedure. The network shown here consists, before prun¬ 
ing, of 800 fibers with Iq = 0.9, the radius and length 
of the cylinder are R = 7 r _1 and L = 10. The filaments 
are fixed at the ends of the cylinder. The effective spring 
constant Y = 1, and k = 10. One may glean from Figs 
3(al),(a2),(bl) and (b2) (which show the network on the 
cylinder, and rolled out onto the (it, v)-plane) in particular, 
there is considerable reorientation in the relaxed network, 
which has an overall energy that that is a factor 2.5 lower 
than the initial network. Several filaments have been plot¬ 
ted in red, to illustrate that initially straight filaments 
are now bent. Also, it is apparent from the increased size 
of the open spaces that - while the overall density has not 
changed - there is considerable stretch and compression 
present in the relaxed network. The most striking change, 
however, is manifested in the orientational distribution 
P(0), which gives the probability of a segment’s angle 
with respect to the long axis being equal to (f>, normalized 
over the interval —7r/2 < < 7 t/ 2 . Where initially, this 

distribution is uniformly flat (P(<(>) = by definition 
for an unstrained mikado network), it acquires a marked 
bimodal structure, with two peaks symmetrically located 
around <j> = 0, see Figs 3(a3) and 3(b3). For lack - to our 
knowledge - of an existing term we shall dub this type of 
order cross-hatched. It is distinct from the more common 
tetratic order, reported for several liquid crystalline sys¬ 
tems DSHD]. Tetratic phases, mostly seen in lyotropic 



FIG. 3: (al) The initial mikado network, wrapped around 
the cylinder. (a2) the same, rolled out onto the ( u , v)-plane. 
(a3) the initial angle distribution P {<j>) is flat, (bl,2,3): same 
images, for the network configuration that minimizes Eq. [3] 
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systems, are composed of two distinct types of nematic 
domains with local directors at right angles to each other. 
A true tetratic phase has fourfold rotational symmetry. 
Though the tetratic phase represents a specific instance of 
cross-hatched order, it possesses a higher symmetry than 
the general cross-hatched phase, which has two superposed, 
nematic phases, with two directors that are generally not 
oriented at right angles to each other and thus has two 
separate twofold rotational symmetries. 

To quantify the extent of the cross-hatched order (de¬ 
noted CH, in the following), we introduce an order pa¬ 
rameter that respects the symmetries of the CH state. 
First, we define the CH-director as the unit vector 
h = A/|A| = ( n u ,n v ), with 

„ r */2 

A = / (ucos |0| + vsin |0|) P(0)d0. (5) 

J-tt/2 

The CH director points along the direction of the positive- 
cA peak in P( 0 ) for sufficiently narrow distributions, and 
is oriented at an angle ( 0 ) = cos _ 1 (n u ) with respect to 
the cylinder’s long axis. We now define the CH-order 
parameter S x to be 

s X=[ cos(4(|0| - (0))) P(0)d 0. ( 6 ) 

J —tt/2 

This order parameter S x is zero in an isotropic network, 
and 1 for any perfectly cross-hatched network - i.e., a net¬ 
work where P(0) = (1/2)[<S(0 —0o) + (S(0 + 0o)], regardless 
of the value of 0o- Similar to its conventional nematic 
counterpart (which is zero for a CH phase), S x quantifies 
the extent of the CH order, but does not contain informa¬ 
tion regarding its spatial orientation. Clearly, the initial 
network as depicted in Figs. 3(al) and 3(a2) has S x = 0, 
the relaxed network depicted in Figs 3(bl) and 3(b2) has 
S x = 0.50. Thus, the curvature of the substrate is seen 
to evoke, in this case, significant cross-hatched order. 

To determine the dependence of the curvature-induced 
CH order and orientations on the radius of the cylinder, we 
perform a series of simulations of the network on cylinders 
of increasing radius. The results are collected, for a 
representative choice of network parameters, and averaged 
over an ensemble of 20 initial networks, in Fig. [4] The first 
striking observation is the emergence of strongly bimodal 
orientational distribution at small cylinder radii. This 
is easily seen in P (</>), and likewise reflected in the order 
parameter. For large radii, we find disordered networks 
with S x = 0, and (0) = 7 r/ 4 . For smaller radii, we see 
S x come up, as well as the emergence of average angles 
(0) smaller than 7 t/4 (note, that a given value of (0) 
corresponds to two peaks in P, on either side of 0 = 0). 
We see that the network is longer in an ordered state as 
Y increases. 

Conclusions. In this Letter, we have studied the spa¬ 
tial organization of two-dimensional filament networks, 


P(<» P(<M P(<» 



FIG. 4: Orientational distribution functions in energy- 
minimized networks on cylinders of radius (a) R = 1/n = Ro, 
(b) R = 6.1-Ro, (c) R = 41Y.Ro. The double peak is seen to 
dissappear at larger radii, as expected. In (d) and (e) the 
cross-hatched order parameter Sx and (0), which is the orien¬ 
tation of the peak for highly ordered networks, is shown. The 
black, blue and green are the data from simulations with 1, 
10 and 100 for Y, respectively and k = 1. 


confined to the surface of a cylinder. We show that the 
curvature acts as a field along the axial direction. A 
simple ring model reveals that this field brings about a 
bifurcation in the optimal orientation of a polymer that 
encircles the cylinder: below a critical radius R *, two de¬ 
generate, slanted ground states emerge whereas for radii 
larger than R *, only the orientation perpendicular to the 
cylinder axis is stable. When instead 2D networks of 
polymers are wrapped around cylinders, this removes the 
sharp bifurcation structure of the ring model, but its gen¬ 
eral characteristics persist: ordered, cross-hatched states 
appear at smaller radii, where the effects of curvature 
along the substrate are largest. The circumferentially 
ordered state is prohibited by the geometry and the fixed 
boundary conditions, and is supplanted by a disordered 
state in which all angles 0 are equally likely. The latter 
is understandable, as at infinite radius the effects of cur¬ 
vature drop out. The effects we observe resemble those 
encountered in various biological settings e.g. collagen 
networks in the growing notochord and in arteries. While 
certainly not the only physical or biochemical principle 
in operation here, we hypothesize that the geometric, 
curvature-induced cross-hatched order we report here 
may explain, in part, the emergence of similar symme¬ 
tries in filamentous biomaterials. In addition, the effects 
we report here may be used to induce order in synthetic 
or biomimetic systems: The instructive effects of curva¬ 
ture could provide a novel and noninvasive route towards 
preparing substrates to predetermined orientational spec¬ 
ifications, and we are currently working to extend the 
framework presented here to arbitrary curved geometries. 
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